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Understanding glassy dynamics is one of the greatest 
challenges in theoretical physics. Many of the recent de- 
velopments in this field are based on the analysis of mean- 
field models 0, Q. The latter are defined by Hamilto- 
nians with long-range interactions and seem to capture 
many properties of real systems. The rather accurate 
comparison to numerical simulations d and experi- 
ments 0, 0, 0| supports the claim that the mechanism 
in these models is similar to the one responsible for the 
glass transition and the glassy dynamics in real materi- 
als. 

In this paper we will consider the so-called Hamiltonian 
Mean Field (HMF) model, a system of planar rotators 
originally introduced in Ref. |8j . This model has been in- 
tensively studied in the last years for its extreme richness 
and flexibility in exploring the connections between dy- 
namics and thermodynamics in long-range many-body 
systems. In fact, on one hand the model has an exact 
equilibrium solution, on the other hand, because of the 
presence of a kinetic energy term in the Hamiltonian, 
the dynamics can be studied by means of molecular dy- 
namics simulations [1 H 0, [Tl| . From these investiga- 
tions, many new interesting features have emerged which 
are common to other systems with long-range interac- 
tions 0, 0, 0| . One of the most intriguing character- 
istics of the dynamics is the existence of quasi-stationary 
states (QSS), i.e. metastable dynamically-created states, 
whose lifetime diverges with the system size N [15|. In 
such states, that spontaneously appear in the numeri- 
cal simulations when the system starts from strong off- 
equilibrium initial conditions, many anomalies have been 
observed, such as anomalous diffusion 0] , non-Gaussian 
velocity distributions , vanishing Lyapunov exponents 
Ra, weak ergodicity breaking, hierarchical structures 
[lq. slow-decaying correlations and aging [l7l ITsl Il9ll20| . 
These features have suggested a possible application of 
Tsallis generalized thermodynamics HI El Hi S El 
but also an interesting link with glassy dynamics. 
We have shown in previous papers that, in the QSS 
regime, the HMF system behaves very similary to a 
spin-glass system [2! [2(j. Actually, by means of 
a new order parameter called polarization and inspired 



by the Edwards- Anderson spin-glass order parameter 

Hi mm ill, 

it has been possible to characterize the 
dynamically generated QSS as a sort of glassy phase of 
the HMF model, despite the fact that neither quenched 
disorder nor frustration are present in the interactions. 
In this paper, by means of a Replica-Symmetry analysis 
performed on an appropriate effective Hamiltonian, we 
will show that it is possible to find out a self-consistent 
equation for a spin-glass order parameter describing, in 
the thermodynamic limit, the quenched dynamics ob- 
served in the QSS regime. We will also show that the 
solutions of this equation reproduce well the microcanon- 
ical simulations results for the polarization in the energy 
region where the dynamical anomalies are more evident, 
thus strongly suggesting the identification of the two or- 
der parameters and confirming the interpretation of the 
limiting QSS regime as a glassy phase. 



I. SPIN-GLASS MODELS 

In the last three decades spin-glasses have attracted the 
attention of experimentalists and theoreticians as glassy 
prototypical systems showing frustration and quenched 
disorder [H El HI HI Hi- Shortly, spin-glasses 
are systems with localized electronic magnetic moments 
whose interactions are characterized by quenched ran- 
domness: a given pair of spins have a roughly equal a 
priori probability of having a ferromagnetic or an anti- 
ferromagnetic interaction. The prototype material is a 
dilute magnetic alloy, with a small amount of magnetic 
impurity randomly substituted into the lattice of a non- 
magnetic metallic host. In this situation, the impossibil- 
ity to minimize simultaneously the interaction energies of 
all the couple of spins leads to a frustration which deter- 
mines a very complex energetic landscape in phase space. 
The latter appears as consisting of large valleys separated 
by high activation energies. Each valley contains many 
local minima in which the system, at low temperature, 
can remain trapped for a very long time. This time grows 
exponentially with the height of the energy barriers, thus 
the system shows very slow relaxation, strong memory ef- 
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fects and aging. 

The modern theory of spin glasses [53 began in 1975 with 
the work of Edwards and Anderson (E A) [23 , who pro- 
posed that the essential physics of spin glasses lays not 
in the details of their microscopic interactions but rather 
in the competition between quenched ferromagnetic and 
antiferromagnetic interactions (i.e. in the frustration). 
Thus they proposed a short range simplified model for 
spin-glasses, in which one represents the magnetic impu- 
rities with Ising spins Sj = ±1 placed on the vertices of 
a three dimensional cubic lattice. The random nature of 
the interactions are mimicked with first neighbors ran- 
dom interactions between the spins taken from a Gaus- 
sian probability distribution with zero mean and variance 
Jfj = J 2 /(2z) where z is the connectivity of the lattice. 
The Hamiltonian (in the absence of an external magnetic 
field) is 

H AS] = ~E JijSiSj , (1) 

(ij) 

where the vector S encodes the full set of spins in the 
sample S — (si, S2, . . . , sjv) and (ij) represents nearest 
neighbors on the lattice. Shortly after the appearance of 
the EA model, an infinite-ranged version was proposed 
by Sherrington and Kirkpatrick (SK) |2SL l30j . For a sys- 
tem of N Ising spins, and in zero external field, the SK 
Hamiltonian is 

Hj[S]=-^=J2J ij s i s j , (2) 

{id) 

with a Gaussian distribution of interactions. Notice that 
in Eq[21 the sum runs over all pair of spins and the fac- 
tor l/y/~N allows to consider the thermodynamic limit 
for the free e nerg y and for other thermodynamic quan- 
tities. In ref.[3Jj| SK showed that their model has an 
equilibrium phase transition to a spin-glass phase below 
the temperature T c = 1 and for an opportune choice of 
the parameters Jo and J, respectively mean and standard 
deviation of the Gaussian distribution of interactions. 
A mean field theory, employing the Onsager reaction field 
term, was proposed two years later by Thouless, Ander- 
son, and Palmer j^], which indicated that there might 
be many low-temperature solutions corresponding to dif- 
ferent spin-glass phases. But the correct solution for the 
low-temperature phase of the SK model, due to Parisi 
[23, employed a novel ansatz and required several more 
years before a physical interpretation could be worked 
out |4{j ■ The picture that finally arose was that of a sys- 
tem with an extraordinary new kind of symmetry break- 
ing, known today as " Replica Symmetry Breaking" , or 
RSB, after the mathematical procedures used to derive 
it. The essential idea is that the low-temperature phase 
consists not of a single spin-reversed pair of states, but 
rather of "infinitely many pure thermodynamic states" 
|40| , not related by any simple symmetry transformation. 
Within the original mean-field framework of the SK 



model it is possible to observe three different phases, 
namely, paramagnetic (PA) , ferromagnetic (FE) and spin 
glass (SG) phase, depending on the temperature and the 
parameters of the Gaussian distribution of the interac- 
tions [2^ . Each phase is characterized by a different mi- 
croscopic behavior and a different kind of orientation or- 
der, giving rise to a characteristic value of the usual mean 
field order parameter, i.e. the magnetization m(T). Thus 
it is clear that the magnetization, calculated at one in- 
stant of time, vanishes in the SG phase just like in the 
PA one. Therefore, in order to discriminate between spin 
glass disorder and paramagnetism, one needs an addi- 
tional order parameter. Such a parameter was originally 
proposed in refs. [13, H3- It is called 'Edwards- Anderson 
(EA) order parameter' and takes into account the tem- 
poral evolution of each spin. In this way the latter is able 
to measure the degree of freezing of the system and to 
distinguish between the PA phase, where it vanishes to- 
gether with the magnetization, and the SG phase, where 
it remains finite. Actually, in its original formulation, the 
EA order parameter results to be only an approximation 
(the so-called Replica Symmetry Approximation) of the 
true spin-glass order parameter proposed by Parisi |39| . 
that has to be defined in the 'replic a sp ace' and results 
to be a whole function (see also refs. [33,133 for more de- 
tails). However, we will show that also considering the 
approximated SK point of view it is possible to shed new 
light on the already suggested link between glassy dy- 
namics and the anomalous out-of-equilibrium behavior 
of the Quasi Stationary States of the HMF model. 



II. DYNAMICAL FRUSTRATION AND 
POLARIZATION IN THE HMF MODEL 

The HMF model describes a system of A*" fully-coupled 
classical inertial XY spins with unitary module 

l?i= (cos Oi, sin 9i) i=l,...,N . (3) 

For a better visualization these spins can also be imag- 
ined as particles rotating on a unit circle. The equations 
of motion follow from the Hamiltonian 

N 2 N 

^E^ + ^Et 1 - - . ( 4 ) 

i—l ij — 1 

where Oi (0 < Oi < 2tt) is the angle and pi the conjugate 
variable representing the rotational velocity of spin i. 
The equilibrium solution of the model in the canonical 
ensemble predicts a second order phase transition from 
a high temperature paramagnetic (PA) phase to a low 
temperature ferromagnetic (FE) one 0, 0, 0, El ■ The 
critical temperature is T c = 0.5 and corresponds to a 
critical energy per particle U c — E c /N = 0.75. The 
order parameter of this phase transition is the modulus 
of the average magnetization per spin defined as: M = 
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FIG. 1: In this figure we plot - for N=100000 - the molec- 
ular dynamics numerical results corresponding to QSS (open 
squares), compared with both the canonical caloric curve, 
plotted as solid line, and the zero magnetization (or mini- 
mum temperature) line, reported as dot-dashed. The tem- 
perature is calculated by means of the average kinetic energy 
K,i.e. T = 2K/N. The QSS follows the latter line only up 
to U = 0.68, becoming unstable below this limiting value 

E3DJ 



wl£i=i ~s~i | ■ Above T c , in the PA phase, the spins 
point in different directions and M ~ 0. Below T c , in 
the FE phase, all the spins are aligned (the rotators are 
trapped in a single cluster) and M ^ 0. 



However, as already pointed out in the introduction, 
molecular dynamics simulations, at fixed energy and for 
out-of-equilibrium initial conditions, show the presence 
of long-living quasi-stationary states (QSS) in the sub- 
critical energy density region 0.5 <U <U C (HHQ1). 
In the thermodynamic limit - when they become sta- 
tionary and reach a limiting temperature value Tqss{U) 

- the QSS can be usefully divided in two subset of states, 
depending on their magnetization: 

- M = QSS (for 0.68 < U < U c ), i.e. anomalous 
states lieing, for N — > oo, over the extension of the high 
temperature branch of the caloric curve (see Fig^l and 
showing a vanishing magnetization (M = 0). It should be 
stressed that, because of the energy density-temperature 
relationship (@) : 

U=^+ 1 -{l-M*) , (5) 

the zero magnetization constraint forces the system to 
follow the minimum kinetic energy (maximum potential 
energy) line, i.e. the minimum temperature line (dot- 
dashed in Fig|TJ; 

■M^0 QSS (for 0.5 < U < 0.68), i.e. anomalous states 
with a macroscopic magnetization (M ^ 0), which for 



N — > oo tend to rejoin the canonical equilibrium curve. 
It has also been shown (H that the M=0 QSS fam- 
ily is unstable below the limiting energy density value 
U ~ 0.68, corresponding to a temperature Tqss = 0.36, 
around which the anomalies are more evident (this is also 
the reason why, historically, the QSS analysis focused on 
the value U = 0.69). Below such a value the QSS cannot 
follow the zero magnetization (minimum temperature) 
line and the dynamical anomalies start to decrease until 
they completely disappear below U = 0.5. 
Actually, it was just the vanishing magnetization of the 
M=0 QSS and the discovery of aging 01 and of dy- 
namical frustration [l^, i.e. the formation of clusters of 
particles with power-law size distributions that sug- 
gested the interpretation of such a regime in terms of a 
sort of spin-glass phase characterized by an EA-like order 
parameter |25j|. 

During the QSS regime, the mean-field interaction seems 
to be broken: in fact, the clusters that compete in trap- 
ping the particles on the unitary circle, generate a dy- 
namically frustrated scenario, in wich each particle never 
feels all the other particles, but only a restricted number 
of them. This feature seems to indicate, for finite sizes 
of the system, a corrugated potential landscape for the 
single particle. This effect is very sensitive to the initial 
conditions |lt| . In the thermodynamic limit (when the 
QSS become stationary), the force between the particles 
- which depends on M - vanishes and the dynamics is 
quenched (apart from the global motion imposed by the 
conservation of total momentum) . 

Therefore, inspired by the physical meaning of the EA 
order parameter, we proposed [2^ a new order param- 
eter for the HMF model, the so-called polarization, in 
order to measure the degree of freezing of the particles 
(rotators) in the QSS regime and, thus, to characterize 
in a quantitative way the emerging glassy-like dynamical 
frustration. 

Polarization is defined as the following spatial average 
1 N 

P = n £ 1 <Tt> 1 ' (6) 
i=i 

where 

i rto+r 

<sl>= - / Ti (t)dt i = 1,...,N (7) 

is the elementary polarization, defined as the temporal 
average, integrated over an opportune time interval r, of 
the successive positions of each rotator. 
For the typical energy density value U = 0.69 (cor- 
responding to a limiting temperature Tqss — 0.38), 
we showed that, while the magnetization correctly van- 
ishes in the thermodynamic limit, the polarization re- 
mains approximatively constant and different from zero, 
thus quantifying the freezing of the rotators in the QSS 
regime. On the other hand, as shown in ref. [25ll26l|. for 
U > U c the polarization coincides with the magnetiza- 
tion and goes to zero. 
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In other words, the polarization seems to play here the 
same role played by the Edwards- Anderson [2^, |3(| or- 
der parameter q^A in the SK model, thus characterizing 
the anomalous QSS regime, which becomes stable in the 
thermodynamic limit, as a sort of Spin- glass phase for 
the HMF model. 

In FigEl extending the results reported in ref. p^LEfil] 
for U — 0.69, we plot, as open circles, the values of the 
polarization order parameter for the range 0.68 < U < 
U C (U C = 0.75), where the M=0 QSS family result to be 
stable. These points, reported for convenience as a func- 
tion of the corresponding N — > oo limiting temperature 
Tqss (0.36 < Tqss < T c = 0.5) were obtained by means 
of the usual microcanonical molecular dynamics simula- 
tions for iV = 1000. The standard "water-bag" initial 
conditions with initial magnetization M — 1 (all the an- 
gles equal and velocities uniformly distributed accord- 
ing to the available kinetic energy) were considered. As 
usual, the integration time r (r = 2000) has been cho- 
sen in order to stay inside the QSS plateaux for every 
N [lit I25L l2fil ] . An average over 25 events was also per- 
formed. The error bars represent the fluctuations of the 
elementary polarization over the configuration of the N 
rotators. As previously seen for U = 0.69 0, increasing 
the size of the system the average polarization remains 
almost constant inside this error along the QSS plateaux. 
Please note that for 0.3 < U < U c the system is in the so- 
called translational regime and its center of mass drifts 
0, i.e. the system is always moving with a global re- 
sulting motion that can be expressed by the phase of the 
average magnetization. Then, as done in the previous 
calculations [25t [2(| , in order to compute the elementary 
polarization, one has to subtract this phase from the spin 
angles. As expected, the polarization correctly decreases 
approaching the phase transition, after which, in the ho- 
mogeneous phase (not shown), it vanishes - together with 
the magnetization - in the thermodynamic limit. 
In the following sections we give further support to the 
claim that the glassy-like behavior characterizing the 
QSS regime corresponds to a 2-vector SK Spin-Glass 
phase by means of a Replica Method formalism applied 
to an appropriately chosen Hamiltonian. 



III. THE REPLICA METHOD 

Suppose |33j that we have a system characterized by an 
Hamiltonian Hj[~s] depending on the configuration [~s] 
of the N spins and on some quenched variables J's chang- 
ing on a time scale infinitely larger than the ~s*'s. If we 
also suppose that these control variables are distributed 
according to a given probability distribution P[J], for 
each choice of the J's one can calculate the partition func- 
tion 

Zj = Y / ^P(-PHj[t}) , (8) 
{«} 
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U=0.68 U=0.73 U =0.75 
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FIG. 2: Microcanonical simulations for the polarization, 
Eq.(6), performed in the QSS regime of the HMF model 
(open circles), are compared with the solution of the self- 
consistent equation for the spin-glass order parameter q (full 
line) Eg. 13011 . Both these quantities are plotted versus the 
limiting temperature Tqss- In the inset, a plot of the caloric 
curve and the QSS points for the HMF model is reported for 
comparison. The spin-glass phase region has been framed. 

and the free energy density 

The point is that averaging the free energy density over 
the distribution P[J], i.e. 

77 = E p [ J ]/^ ■ ( 10 ) 
j 

results to be a task not simple at all. Thus in refs p^Esj 
the so called Replica Method was proposed. The latter is 
a trick to simplify the calculation of Eq. and consists 
in computing the average of the free energy density by 
some analytic continuation procedure from the average 
of the partition function of n uncoupled replicas of the 
initial system . In fact, using the identity 

lnZ, 7 =lim i^-ll , 

n— >0 n 

togheter with Eq.@, Eq. i|10|l can be rewritten as 

j 

= -/3- 1 Hm lim-^£P[J]{(^r-l} . (11) 
J 

Finally, if we define 

(z7r = j2p[J]{Zj} n (12) 
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and we use the normalization condition Y] T P\J] = 1, 
Eci. (|ll[) becomes 

f = -P~ 1 Mm lim-UpjF-l} (13) 

AT— »oc n— >0 niv 

Now, denoting with a the replica index (a = 1, ...n, with 
integer n) , we can write the partition function in Eq. ({12(1 
as the partition function of n non-interacting replicas of 
the same system (for the same set of J's) 

n 

W = E E " E exp{-^/3i/. 7 [^]} 

{ S 1}{ S 2 } {s n } a = l 



= Tr cxp{ - £ PHj{l? a }} , (14) 

a=l 

where the trace Tr in the last expression syntcthizes the 
sums over the spins in all the replicas. Further averaging 
this quantity over the distribution P[J], the calculation 
of the averaged free energy density follows straightfor- 
ward from Eq. (|13l) . 

Our idea is to describe the glassy dynamics of the M=0 
QSS regime of the HMF model by studying the equilib- 
rium properties of an infinite range XY spin-glass effec- 
tive Hamiltonian 

Hj[~s] = -gE ■ ~$ ' ( 15 ) 

id 

with an opportune choice of the interactions J,j. 
In Ea. (|15|> each spin (rotator) is defined as = 
(cosQi, sinOi), with < 9i < 2n and unitary module, 
while the factor ^ before the summation prevents from 
counting two times the same spin couples. The distribu- 
tion of the quenched variables Jij has to be chosen such 
that: 

• it must take into account the presence of the dy- 
namically created clusters of particles observed in 
the QSS regime, which in turn generate dynamical 
frustration; 

• its first moment J a should be equal to 1/N (being 
in the HMF model = 1/N Vi, j); in such a way, 
for N — ► oo, Eq. fTSj) will reduce to the potential 
term of the HMF model. 

Thus, without loss of generality, we can choose a Gaus- 
sian distribution: 

p(J l3 ) = [^/V]- 1 e X p - {J Y J2 J0? (16) 

where the first two moments, the mean J a and the vari- 
ance J 2 , will be set equal to 1/iV. 

Finally, the inverse temperature /3 = 1/T will be fixed 
by the kinetic term of the HMF model Hamiltonian (be- 
ing T — 2K/N), that we assume to be constant for the 



canonical procedure we want to perform. Of course we 
are implicitly assuming that the stationarity of the M=0 
QSS regime in the thermodynamic limit could enable us 
to use equilibrium thermodynamics tools. 
For fixed J's we expect that the replica method would en- 
able us to find out a self-consistent equation for the spin- 
glass order parameter of the model l|15f) . This equation 
will be obtained using the Replica- Symmetry (RS) ansatz 
in the context of a 2- vector infinite range spin-glass model 
|30| , by imposing the spin-glass extremal constraints dur- 
ing the steepest descent procedure. Our final goal will be 
to compare such a theoretical prediction with the molec- 
ular dynamics results for the polarization of Eq.(jSJ, in a 
range of temperatures corresponding with the N — * oo 
limiting temperatures (Tqss) that characterize the ho- 
mogeneous quasi-stationary states. 

Let us start by applying the replica trick to the free en- 
ergy calculation starting from the effective Hamiltonian 
()15p. with quenched couplings Jij following the distribu- 
tion l(ISjL 

Using some general properties of the characteristic func- 
tion of a statistical distribution, and averaging Ea. l|14l) 
over the chosen distribution of the , we can write the 
free energy density expression (|13|l as a function of only 
the first two moments of the distribution itself 

f = — B^ 1 lim lim 

N^oo n->0 nN 
(ij) a 

a b 

where a, b are replica indexes and i,j are spin indexes 
in each replica. Thus ~Si a = (cosOi a , sinO®) represents 
the ith spin in the ath replica. As seen before, both the 
parameters J and J 2 , respectively mean and variance of 
the J's distribution 116|l . have to be put equal to 1/N. 
Finally, the notation in the sum is equivalent to a 
sum over all N(N — l)/2 distinct pairs of sites, thus the 
factor 1/2 disappears behind the sum over i, j. 
Because of the latter notation, and after some rearrang- 
ment, we can write the following equivalences 

EE^-^EiE^i 2 -^ ( 18 ) 

(ij) a a i 

and 

/] El S{ ' s j El Si ' s j 

(ij) a b 

-EE cos ^ a -^) cos ^ b -^) 

a,b (i,j) 
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= E E \ ' + ■ Tj ab ] (19) 

where two terms of interference between replicas appear 
Sl ab = (cos^ - 0,: b ), S m(^ a - Of)) , 

and 

i\ ab = (cosies + eS), sinies + e b )) . 

The latter term, for a = b, becomes US — 
(cos29S,sin29S). 

After further rearrangment in the summations, we obtain 
the following expression for the free energy density 



1 



f = —0 lim lim iexp 

AT^oo n ->0 nN 



P -(nN~2n>)- n -E 
8 2 



xTrexp [ AvlV^ a|2 



2N 

a i 



8N 



with the positions 



T 



_ \ ^ 7ft ab ~Qab _j_ \ ^ ~Q*ab 

i i 

thus obtaining 

I = Tr f ]J(— ds^du a ) Y[(— dS ab dt ab ) 

J a 77 ab 77 

xexp {-iV[i^|s<f + i^|#f 



\ Ed^l 2 + I*"!*) " & J2s°~S a ' ^ a U a 



a^b 



£j2(s ah ^ ab +? b 3 tob ) ] } 



EiE^ a i 2 + E(iE^ ab i 2 + iE^ a 

i a^b i i 



b\2\ 



i} 



(20) 



In the thermodynamic limit, the first exponential be- 
comes 

r/? 2 / , T ~ 2, n/3. r nN0 2 ^ 
exp [—(nN -2n 2 ) -] w exp [ — ] 



The mean field approximation consists on saying that 
one can deal with independent spins feeling the external 
conjugated fields ■ so that this last relation becomes 

J a ab 

x Tr n [ exp {-\ £ |?f - i £ |«»| a 



and it does not involve the glassy properties of the sys- 
tem. Therefore in the following we will concentrate on 
the term 





/ = Trexp t^EiE^ a 



8N 



EiE^ Q i 2 + E(iE^ ab i 2 + iE^ a 

i a^b i i 



b\2\ 



(21) 



It can be linearized, term by term, with the Hubbard- 
Stratonovich (HS) Gaussian transformation, i.e. 

exp [nx 2 ] (2-kY 1 ! 2 I dyexp [- V — + (2^ 2 xy } , (22) 



4Ed 



grab|2 _|_ l^txb 1 2 \ _|_ fli \ ~ ~c a , ' Z7 a Tr a 



p^^s* S a + ~v?U l 



a^b 



(3 



J2(S ab ^ ab + t ab T ab ) } 



(23) 



a^b 



where from now on the trace is over the n replicas at a 
single site (mean field approximation). 



IV. SPIN-GLASS SELF-CONSISTENT 
EQUATION AND ITS NUMERICAL SOLUTION 

At this point we have to impose the following spin- 
glass constraints [3(j, based on the physical meaning of 
the two order parameters in the M=0 QSS regime (that 
we are considering as an effective spin-glass phase of the 
HMF model): 
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1. the first one refers to the magnetization, that is 
null in the M=0 QSS regime, thus implying 



(24) 



2. the second one refers to the spin- glass order param- 
eter, which of course does not vanish in the spin- 
glass phase thus quantifying the degree of freezing 
of the rotators in the M=0 QSS regime; as usual 
in standard replica method for glassy systems, it is 
chosen as proportional to the module of the overlap 
between two different replicas at a single site 

q n = q ab = 2[3- 1 \r b \ , (25) 

where the |s° h | are considered equals Va, b (Replica- 
Symmetry approximation). The phase of s" , being 
arbitrary, can be set to zero for convenience so that 
its direction is u£. Finally, we set also | t ab \ =0. 

By incorporating the spin-glass constraints with delta 
functions one can easily perform the integral in Ea. (|23[l . 
thus obtaining 



/ = [ Tr n exp( - ]T 



4 ^4 

a^b a^b 



2 , — — u x s 



ab \ lJV 



(26) 

Since the first term is independent on the replicas and 
the trace being a linear application, one has equivalcntly 



a=tb 



r°° r 2 r>n ft 2 , 

-In J rdrexp[----^-] V(/WW2) ] ) . 

(29) 

The latter exponential, for N — > oo, will show a maxi- 
mum at the extremum of its argument, thus we have to 
impose the following (steepest descent) extremal condi- 
tion 

d , n(n - 1) (3q n 2 
dq n [ 2 1 2 ' 



-In / O °rdrexp[-^-^^]V( / 3r v ^72)]=0 . 
JO Z 4 

Performing the derivative we find the desired self- 
consistent equation for the order parameter q n 



q n (l - n) = 1 



J^r 2 dreM~ r -} Io n (/3r^/2) f 



/ °°rdrexp[-^] I n ((3ry/qj2) 
that has to be continued for n — > 0, giving finally 

g = l-JV f\^exp[-i^ 

V q Jo 2 j 



(30) 



exp{ — N[ 



n(n - 1) (3q n 



-lnTr n exp(-^E lI *~S ab )}) • (27) 

a^b 

Reminding that S ab is the replica interference at single 
site, i.e. ~S ab = (cos{8 a - 8 b ), sin(8 a - 8 b )), one has that 
~Ux S ab = cos(8 a — 8 b )~s a ~s b , so we can write 



EuJ ab = C£^ a ) 2 - 

a^b a 



(28) 



Inserting the latter relation in Eq. l|27[) . using 
the Hubbard-Stratonovich transformation (for 

exp[— 4^(Y2 a ~s a ) 2 ]) onc e again and substituting 
the trace with the following multiple integral 

m..,-rii £[...], 

a— 1 



(I 



one finally finds 



t < at M n ~ 1) ,Pq n , 2 
I = exp {-N [ (— ) 



Solving numerically this equation, we can immediately 
calculate the expected value of the spin-glass order pa- 
rameter q as a function of the temperature T = 
and compare it with the polarization order parameter 
p of the HMF model. In particular, since we have to 
consider M~0 QSS with null magnetization, we consider 
the limiting temperature Tqss (i.e. the temperature for 
N — ► oo) and, by varying it in Eq. i|3Ufl . we obtain 
the theoretical curve shown in Fig[21 (full line). This 
curve is compared with the molecular dynamics simu- 
lations for the polarization (open circles) performed in 
the QSS regime of the HMF model for N=1000 and for 
different temperatures in the subcritical region. Some 
corresponding energy density values are also reported for 
convenience. As previously stressed, these results are in- 
dependent on N within the error bars ^(| and this allow 
us to extrapolate them to the thermodynamic limit. 
The theoretical curve predicts, as expected, a phase tran- 
sition at T c = 1/2 and superimposes on the values of the 
polarization obtained in the range 0.36 < Tqss < 0.45, 
that corresponds (being there Mqss — 0) to the energy 
density range 0.68 < U < 0.72 through the zero magneti- 
zation (minimum temperature) line equation U — -j + \ , 
see Eq.(|5}. The simulations points start to disagree with 
the results of the theoretical curve around Tqss ~ 0.45, 
i.e. for U > 0.72: above these values the QSS points 
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are very close to the correspondent equilibrium temper- 
ature values on the caloric curve (as visible in the inset 
of FigEJl and the glassy features of the QSS regime tend 
to disappear. Inverting the previous argument, we could 
also suggest that such a theoretical result allows us to 
better specify the range of energy densities where the 
QSS regime can be considered as a spin-glass phase for 
the HMF model. 



V. CONCLUSIONS 

We have discussed the glassy phase of the HMF model 
by means of an effective spin-glass Hamiltonian and 
shown that the corresponding order parameter coincides 
with the polarization in the energy range where HMF 
exhibits strong dynamical anomalies and glassy dynam- 
ics. Being a long-range Hamiltonian solvable model and 
showing the existence of a glassy-like relaxation dynam- 
ics along quasi-stationary trajectories, the HMF model 
poses new challenging questions on the origin of its glassy 
unexpected behavior. In these respects, the model may 



be a very useful " laboratory" for studying general trends 
that can be later tested numerically and experimentally 
in more realistic systems and real materials. 
These analytical results, obtained within the Replica- 
Symmetry SK framework and in the thermodynamic 
limit, seem to give further support to the interpretation 
of the QSS regime as a real spin-glass phase. This is 
true at least in the limiting temperature region where 
the glassy features are more evident. In this region the 
elementary polarization seems to play the role of the EA 
order parameter in giving a measure of the numerically 
observed quenched dynamics. 

In conclusion, we do hope that this connection between 
the HMF model and Spin-glass systems could bring new 
insight on the several common features. 
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